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Abstract 



We have fully reconstructed decays of both B° and B~ mesons into final 
states containing either D, D* , D**, ip, if)' or Xd mesons. This allows us to 
obtain new results on many physics topics including branching ratios, tests 
of the factorization hypothesis, color suppression, resonant substructure, and 
the B~ — B° mass difference. 
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I. INTRODUCTION 



Since B mesons were first fully reconstructed in f 983 by CLEO [f ] there have been several 
papers by CLEO [2], [3] and ARGUS [4], [5], [6] which reported branching ratios for exclusive 
decay modes of B mesons. We present here new data from the CLEO II detector using a 
high resolution photon detector and a much larger data sample than has been available 
previously. 

We are particularly interested in two-body hadronic B meson decays, which occur through 
the Cabibbo favored 6 — > c transition. In these circumstances the dominant weak decay 
diagram is the spectator diagram, shown in Fig. 1(a). The virtual W~ materializes into 
either a ud or cs pair. This pair becomes one of the final state hadrons while the c quark 
pairs with the spectator anti-quark to form the other hadron. The Hamiltonian [7] , ignoring 
hard gluon corrections, is 

H = ^V cb {[(du) + (sc)}(cb)} (1) 

where (q~iqj) = <?;7,u(l — 75)^'? Gf is the Fermi coupling constant, and V cb is the CKM matrix 
element. 

The spectator diagram is modified by hard gluon exchanges between the initial and final 
quark lines. The effect of these exchanges can be taken into account by use of the renor- 
malization group. These gluons induce an additional term so that the effective Hamiltonian 
is comprised of two pieces, the original one now multiplied by a coefficient Ci(fi) and an 
additional term multiplied by c 2 (/i): 

H eff = ^V cb { Cl (fi) [(du) + (sc)} (cb) + c 2 (fi) \(cu)(db) + (cc)(sb)}} (2) 

where the c 4 - are Wilson coefficients evaluated at the mass scale //. The Wilson coefficients 
can be calculated from QCD; however, the calculation of rates is inherently difficult because 
it is unclear at what scale these coefficients should be evaluated. The usual scale is taken to 
be fi ~ ml. Defining 

c±(aO = c i(/«) ± c 2(/«) (3) 
the leading-log approximation gives [8] 

-67± 



c±Ul)= » (33 -2,,,) 

V <w) / 

where 7_ = — 27+ = 2, and nj is the number of active flavors, five in this case. 

The Hamiltonian in Eq. (2) leads to the "color suppressed" diagram shown in Fig. 1(b), 
which reflects the quark pairings in the term multiplied by the coefficient c 2 (/i). Observa- 
tion of B — > ipX s decays, where X s is a strange meson, gives experimental evidence for the 
existence of this diagram. Further information on the size of the color suppressed contribu- 
tion can be obtained from B° — > D° (or D*°)X° transitions, where X° is a neutral meson 
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containing light quarks. In B~ decays, both types of diagrams are present and can interfere. 
By comparing the rates for B~ and B° decays, the size and the sign of the color suppressed 
term can be extracted. 

Bjorken has suggested [9] that, in analogy to semileptonic decays, two body decays of B 
mesons that occur via the external spectator process can be expressed theoretically as the 
product of two independent hadronic currents, one describing the formation of a charm meson 
and the other the hadronization of the ud (or cs) system from the virtual W~ . Qualitatively, 
he argues that for a B decay with a large energy release the ud pair, which is produced as a 
color singlet, travels fast enough to leave the interaction region without interfering with the 
formation of the second hadron. The assumption that the amplitude can be expressed as 
the product of two hadronic currents is called "factorization" in this paper. Several tests of 
the factorization hypothesis can be made by comparing semileptonic and hadronic B meson 
decays. 

This paper is structured in the following manner: the data sample, detector and recon- 
struction procedures are described in sections II and III. Branching ratios are given for 
B — > Dtt~ and B — > Dp~ modes in section IV. In section V results on branching ratios, 
polarizations and final state substructure for B — > D*tt~ } B — > D* p~ and B — > D*a^ are 
described. Section VI describes a search for D** production in hadronic B decay. This is 
followed by section VII on exclusive B decays to charmonium, and section VIII on a search 
for other color suppressed B decays. A B~ — B° mass difference measurement is described 
in section IX. The interpretation of these results and comparisons to theoretical predic- 
tions are discussed in sections X (factorization tests), XI (spin symmetry tests) and XII 
(determination of the color suppressed amplitude). 

II. DATA SAMPLE AND SELECTION CRITERIA 
A. Data Sample 

The data sample used in this paper was collected with the CLEO II detector at the 
Cornell Electron Storage Ring (CESR). The integrated luminosity is 0.89 fb -1 at the T(4S') 
resonance and 0.41 fb -1 at energies just below BB threshold, henceforth referred to as 
the continuum. It is natural to assume equal production of charged and neutral B } s since 
the difference between their masses is very small (see section IX). Then there are a total 
of 935, 000 ± 10, 000 ± 15, 000 charged and the same number of neutral B mesons in this 
sample. 

B. Detector 

The CLEO II detector [10] is designed to detect both charged and neutral particles 
with excellent resolution and efficiency. The detector consists of a charged particle tracking 
system surrounded by a time-of-flight (TOF) scintillation system and an electromagnetic 
shower detector with 7800 thallium-doped cesium iodide crystals. In the "barrel", defined 
as the region where the angle of the shower with respect to the beam axis lies between 32° 
and 135°, the r.m.s. energy resolution is given by SE/E(%) = 0.35/E ' 75 + 1.9 — 0.1-E (E 
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in GeV). In the endcap region, located between 18° and 36° from the beam axis, the r.m.s 
energy resolution is given by SE/E{%) = 0.26/ E + 2.5 . The tracking system, time-of- flight 
scintillators, and calorimeter are installed inside a f .5 T superconducting solenoidal magnet. 
Immediately outside the magnet are iron and chambers for muon detection. The momentum 
resolution of the tracking system is given by (Sp/p) 2 = (0.0015p) 2 + (0.005) 2 (p in GeV/c). 
Ionization loss information (dE/dx) } provided by the tracking system, is used to identify 
charged particles in this analysis. The track must have a dE / dx measurement that differs 
from that expected for the charged particle hypothesis under consideration by less than 3<r 
(where henceforth a denotes the r.m.s. resolution). 

Muons are identified by a system of drift tubes interleaved with layers of magnet iron. 
Electron identification utilizes the specific ionization of the track in the drift chamber, the 
spatial distribution of the energy in the calorimeter, and the ratio of the cluster energy 
measured in the calorimeter to the track momentum. 

C. Photon Selection 

Photon candidates are selected from showers in the calorimeter barrel that have a mini- 
mum energy of 30 MeV, are not matched to a charged particle track from the drift chamber, 
and have a lateral energy distribution consistent with that expected for photons. In the 
calorimeter endcap the same criteria are applied but the minimum energy requirement is in- 
creased to 50 MeV. A small angular region between 32° and 36° degrees in the barrel-endcap 
overlap region is excluded. Neutral pion candidates are selected from pairs of photons with 
an invariant mass within 2.5<r of the known 7r° mass. These candidates are kinematically 
fitted with a 7r° mass constraint. 

Candidate rj mesons are reconstructed in the ?y — s- 77 mode. They are required to have 
an invariant mass within 30 MeV of the known rj mass (547.5 MeV) [14]. The candidates 
which pass the requirements described above, are kinematically constrained to the rj mass. 
Candidate rj' mesons are reconstructed in the ?y7r + 7r~ channel with rj — > 77. Candidate u 
mesons are reconstructed in the u — > 7r + 7r~7r channel. 

D. Charm meson selection 

We select D°, D + , D* + and D*° mesons based on the following criteria. Candidate 
D° mesons are identified in the decay modes D° — > K~tt + } D° — > K~TT + Tr° and D° — > 
K~ 7T + 7r + 7r~ . Candidate D + mesons are selected using the D + — > K~tt + tt + mode. The 
decay modes, branching ratios and r.m.s. mass resolutions, cr mD , are listed in Table I. We 
use the CLEO [11] absolute branching ratio for D° — > K~tt + decays [13], and the Mark III 
value for D + — > K~tt + tt + [12]. We use the Particle Data Group (PDG) values [14] for the 
ratios B(D° — > K~ 7r + 7r ) / B(D° — > K~tt + ) (where henceforth B denotes the branching ratio) 
and B(D° -> K-Tr + Tr+Tr-)/B(D° -> K~tt + ). 

Charged D* candidates are found using the decay D* + — > tt + D° } while neutral D* candi- 
dates are found using the decay D*° — > tt°D . Other D* decay modes are not used because 
they have much poorer signal to background ratios. CLEO branching ratios (Table II) are 
used for D* decays [16]. We form D* + and D*° candidates by selecting D° candidates whose 
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mass is within 2.5<r of the known D° mass. Then we require that the D*-D° mass difference 
be within 2.5a of the measured values [14,17]. 

E. Charmonium Meson Selection 

We reconstruct the charmonium states ip, if) 1 and xa, where ip mesons are selected by 
their decay into pairs of identified leptons (e + e~ or /i + \i~ ). We use the Marklll value [18] 
B(t/j — > = (5.91 ± 0.25)% for the ip to dilepton branching ratio. The kinematics of 

B decay at the T(4S') imply that each lepton in a ip candidate has a momentum between 
0.8 GeV/c and 2.8 GeV/c, with one lepton always having momentum greater than 1.5 
GeV/c. For the clean modes B~ — > ipK~ and B° — > ipKg, we obtain good efficiency in 
the dimuon channel by requiring only one identified muon that penetrates through three 
interaction lengths. In the dielectron channel, one of the electrons must satisfy a loose 
electron probability requirement. For modes other than B~ — > ipK~ and B° — > ipKg, both 
electrons must be identified, or one muon is required to penetrate five interaction lengths 
and the partner muon is required to penetrate three interaction lengths. 

Final state radiation is included in the Monte Carlo simulation of ip meson decays. For ?/>'s 
in the dielectron final states we employ an asymmetric mass cut: —150 < m(e + e~) — m(ip) < 
45 MeV in order to reduce the efficiency loss from this source. For the dimuon final state, 
we require —45 < m( / u + / u~) — m(ip) < 45 MeV since final state radiation is less significant in 
this case. (The ip mass resolution would be 15 MeV (r.m.s.) in the absence of radiation). In 
these mass windows the efficiency for detecting ip mesons in the dielectron and the dimuon 
final states are 48.1% and 67.8% for the looser cuts, and are 45.7% and 42.4% when both 
leptons are identified. 

The decay modes ip — > e + e~, ip' — > fi + \i~ , and ip' — > 7r + 7r~?/>, are used to select ip' 
candidates. The recontruction of the leptonic decays follows the procedure outlined for ip 
mesons. Pion candidates for the third decay mode are required to have dE/dx measurements 
consistent with the pion hypothesis. In addition, tracks that have been identified as part of 
a Kg decay are rejected. It has been shown that the tt + tt~ invariant mass spectrum from ip' 
decays favor larger values relative to that expected from phase space [18]. We require the 
7r + 7r~ invariant mass to be between 0.45 and 0.58 GeV [19]. For ip' mesons reconstructed 
through the decay ip' — > 7r + 7r~?/> we require the ip — ip mass difference, 8m = m^i — m^,, to 
be between 0.568 and 0.599 GeV. Xa mesons are reconstructed by their decay into a photon 
and a ip meson. We require the photon be in the good portion of the barrel calorimeter 
(| cos#| < 0.71). If the photon candidate forms an invariant mass within —5 to 3 standard 
deviations of the known 7r° mass when combined with any other photon in the event, it is 
rejected. 

III. B MESON RECONSTRUCTION PROCEDURES 

A. Candidate Selection 

After selecting D } D* or charmonium candidates we combine them with one or more 
additional hadrons to form B candidates. The measured sum of charged and neutral energies, 
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Emea, s of correctly reconstructed B mesons produced at the T(4S') must equal the beam 
energy, -Ebeam, to within the experimental resolution. Depending on the B decay mode, a&E, 
the r.m.s. resolution on the energy difference AE = -Ebeam — E meas varies from 8 to 46 
MeV. The modes considered and the corresponding a&E values are given in Tables III— VII 
and Tables IX-X. We divide the B candidates into a signal sample where AE is consistent 
with zero within 2.5 <r, and a "sideband" sample consisting of two intervals, one with AE 
positive and the other negative, both 2.5a wide and at least 3a away from AE = 0. For 
decay modes with large a&E, we restrict the sideband width so that the maximum value of 
AE is less than one pion mass. This avoids contamination from the B decay mode with an 
additional pion. These AE sidebands are used to study the background shape. 

For B decay modes with a fast p~ the energy resolution depends on the momenta of the 
pions from the p~ decay. The momenta of the charged and neutral pions are correlated; a fast 
7r~ accompanies a slow 7r° and vice- versa. This correlation is most conveniently formulated 
as a function of the helicity angle P , the angle in the p~ rest frame between the direction 
of the 7r~ and the p~ direction in the lab frame. When cos P = +1, the resolution in 
the energy measurement is dominated by the momentum resolution on the fast tt~ . In 
contrast, when cos P = —I, the largest contribution to the energy resolution comes from 
the calorimeter energy resolution on the fast 7r°. Typically a&E varies linearly between 20 
MeV at cos P = —I and 40 MeV at cos P = 1. The energy resolution from a Monte 
Carlo simulation for one such mode (B~ — > D°p~) is shown in Fig. 2 as a function of the 
p~ helicity angle. The energy difference resolutions for modes containing a p~ are given in 
Tables III VI. 

In addition to the above selection criteria, events are required to satisfy R 2 < 0.5 
where R 2 is the ratio of the second Fox- Wolfram moment to the zeroth moment determined 
using charged tracks and unmatched neutral showers [20]. A sphericity angle cut is applied 
to further reduce continuum background. The sphericity angle S is the angle between 
the sphericity axis of the particles which form the B candidate and the sphericity axis of 
the other particles in the event [21]. For a jet-like continuum event, the absolute value of 
this angle is small; while for a BB event, the two axes are almost uncorrelated. Requiring 
|cos0 s | < 0.7 typically removes about 80% of the continuum background, while retaining 
70% of the B decays. The sphericity cut used here depends on the number of pions which 
accompany the D or D* meson. For final states with a D* and a single (2, 3) pion(s) we 
require |cos0 s | < 0.9 (0.8, 0.7). For all modes which contain a D and a single (2) pion(s) 
in the final state, we demand that |cos0 s | < 0.8. In modes with ip mesons, we maximize 
the efficiency by applying no sphericity angle cut. 

To determine the signal yield and display the data we form the beam constrained mass 



where pi is the momentum of the z-th daughter of the B candidate. The resolution in this 
variable is about 2.7 MeV [22] and is about a factor of ten better than the resolution in 
invariant mass. The width is dominated by the CESR beam energy spread rather than by 
detector resolution. 




(5) 



For a specific B decay chain, such as B 



D tt- } D° 



K 7T + 7r°, we allow only one 
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candidate per event to appear in the Mb distribution. If there are multiple candidates with 
Mb > 5.2 GeV, the entry with the smallest absolute value of AE is selected. 

B. Background Studies 

In order to extract the number of signal events it is crucial to understand the shape of 
the background in the Mb distributions. There are two contributions to this background, 
continuum and other BB decays. The fraction of background from continuum events varies 
between about 58% and 91% depending on the B decay mode [23]. 

We expect that the Mb distribution from the AE sidebands will give a good represen- 
tation of the background shape. To verify this, a Monte Carlo simulation of BB events was 
used to show that the AE sidebands can be used to accurately model the shape of the BB 
background under the signal in the beam constrained mass distributions (see Fig. 3). In 
continuum data, the AE sidebands also model the shape of the background in the signal 
region. The sum of the BB Monte Carlo and continuum data agrees in shape with the AE 
sidebands in data (see Fig. 4). Therefore, the AE sidebands can be used to model the shape 
of the background under the signal in data. The Mb distributions for AE sidebands in data 
for several modes are shown in Fig. 5. All of these can be fitted with a linear background 
below M#=5.282 GeV, and a smooth kinematical cutoff at the endpoint, which we choose 
to be parabolic. The distributions of Mb for wrong-sign combinations (e.g. D°tt + ) } wrong- 
charge combinations (e.g. D + tt + ) } and continuum data can also be adequately fitted with 
this functional form (henceforth referred to as the CLEO background shape). To determine 
the number of signal events from the Mb spectrum in the AE interval centered on zero, we 
use the background function as determined from the sidebands and a Gaussian signal with 
a fixed width of 2.64 MeV. 

C. Efficiency Studies 

In order to extract branching ratios, detection efficiencies are determined from a Monte 
Carlo simulation of the CLEO II detector. The accuracy of the simulation is checked in 
several ways. We select radiative Bhabha events (e + e~ — > e + e + 7) using only calorimeter 
information and then embed the tracks into hadronic events. We find that the efficiency 
for the detection of charged tracks above 225 MeV/c is correct to better than 2%. The 
Monte Carlo simulation of charged tracks with transverse momenta below 225 MeV/c is 
more complicated since these tracks do not traverse the entire drift chamber. The accuracy 
of the simulation is verified using the D* decay angle distribution of inclusive D* + — > D°tt + , 
D° — > K~tt + decays which must be symmetric after efficiency correction. The simulation of 
low px tracks agrees with the Monte Carlo simulation for 100 < p < 225 MeV/c. However, 
the efficiency for tracks in this momentum range is known to only ±5%. 

The accuracy of the photon detection efficiency can be verified by comparing the ratio 
of branching ratios of 77 — > 7r 7r 7r and ?y — s- 77 to the average ratio given by the PDG 
[14]. This test indicates that the single photon efficiency is modelled to better than ±2.5%. 
Other checks of the tt + and 7r° detection efficiency are performed by comparing the yield of 
fully reconstructed D° — > K~TT + Tr° decays with the yield of partially reconstructed D° — > 
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K 7T°(7r + ) where the tt + is not detected. Additional consistency checks have been performed 
by comparing inclusive D* + and D*° cross sections in the continuum, and by comparing the 
ratios B(D° -> A"-7r + 7r-7r + ) /B(D° -> K~tt + ) and B(rj -> 7r-7r + 7r )/B(?7 -> 77) to the 
values in the PDG compilation [14]. 

IV. BRANCHING RATIOS FOR Dir~ AND Dp~ FINAL STATES 

We reconstruct the decay modes B° — > D + tt~ } B° — > D + p~ , B~ — > D°tt~ } and B~ — > 
D° p~ following the procedures described in sections II and III. There is an additional 
complication for the analysis of the B — > Dp~ modes. Events which are consistent with 
the decay chain B — > D*tt~ } D* — > Dtt° have the same final state particles and thus form 
a potential background. We eliminate this background by discarding events for which the 
D* — D mass difference is consistent with the D* hypothesis. This veto does not reduce 
the efficiency for B — > Dp~ . A Monte Carlo simulation of BB decays shows a broad 
enhancement in the signal region for B~ — > D°tt~ and B~ — > D° p~ (see Fig. 3). This 
enhancement contains contributions from B — > D*°X } D*° — > D°7 transitions which can be 
modeled with the CLEO background shape. 

To select B — > Dp~ channels we impose additional requirements on the 7r~7r invariant 
mass and decay angle. Specifically, we require that |m(7r~7r ) — 770| < 150 MeV/c 2 . Since the 
decay B — > Dp~ is fully longitudinally polarized (helicity zero due to angular momentum 
conservation), a cut on the p helicity angle is imposed (|cos0 p | > 0.4) [24]. The beam 
constrained mass distributions for B — > Dtt~ and B — > Dp~ are shown in Fig. 9. Fig. 6 
shows the p helicity angle distributions (with the cut on the helicity angle removed) for 
B° — > D + p~ and for B~ — > D° p~ after B mass sideband subtraction. After efficiency 
correction, these distributions are given by the functional form: 



where T^/T is the fraction of longitudinal polarization. The fit gives T^/T = 1.07 ± 0.05 for 
B~ — > D°p~ and T^/T = 0.92 ± 0.07 for B° — > D + p~ . These results are consistent with full 
polarization as expected and thus provide a consistency check of the background subtraction 
and efficiency correction. Monte Carlo simulation shows that most of the BB backgrounds in 
B — » Dp~ decays are due to combinations with an incorrectly reconstructed low momentum 
7T°. Therefore a fit to the beam constrained mass distribution with cos P < —0.4 is also 
performed as a consistency check of the analysis [25]. These results agree with the branching 
ratios obtained using the full range of helicity angle. 

The 7T~7r° invariant mass distribution for the B signal region (±6.5 MeV of the nominal 
B mass), is shown in Figs. 7 and 8 after B sideband subtraction. Fitting this distribution 
to the sum of a Breit Wigner and a parameterization of non-resonant B — > Dtt~tt° decay 
[26] we find that fewer than 2.5% (at 90% C.L.) of the events in the B mass peak arise from 
non-resonant decays, after applying the helicity angle cut and restricting the 7r~7r mass to 
lie in the rho mass region. The observed Dtt~tt° events are consistent with B — > Dp~ and 
any non-resonant contribution can be neglected. 
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The resulting branching ratios lor B — > Dtt and B — >■ Dp are given in Tables III and 

IV. 

Two systematic errors are quoted on the branching ratios. The first includes contributions 
Irom background shape (~ 5%), Monte Carlo statistics (2 — 4%), and the uncertainty in the 
modeling of the tracking and 7r° detection efficiencies (which depend on the multiplicity of 
the decay mode as described earlier) and the relative D° branching fractions. The second 
systematic error contains the errors from the D + — > K~tt + tt + (±14%) and D° — > K~tt + 
absolute branching ratios (±2.7%). 



V. MEASUREMENTS OF D*(Ntt)- FINAL STATES 
A. Branching Ratios 

We now consider final states containing a D* meson and one, two or three pions. These 
include the B — > D*tt~ , B — > D* p~ , and B — > D*a~[ decay channels. A cut on the D* 
helicity angle, |cos#d*| > 0.4, is made for B — > D*tt but not for D* p and D*a\. The 
beam constrained mass distributions for the B° — > D* + tt~ and B~ — > D*°tt~ are shown in 
Fig. 10. Our results for the decays B° — > D* + tt~ and B~ — > D*°tt~ are listed in Tables 
V and VI. The first error quoted on the branching ratios is statistical, followed by two 
systematic errors. The first systematic error contains contributions from the uncertainties 
in the efficiency for charged track finding, the uncertainty in photon detection, variations in 
event yield from changes in background shape, Monte Carlo statistics and the relative D° 
branching fractions. The second systematic error contains the errors on the D° — > K~tt + 
and D* — > D°ir branching ratios. 

Fig. 10 shows the beam constrained mass distributions for the B° — > D* + p~ and B~ — > 
D*°p~. To study the resonant substructure in B° — > D* + tt~tt° the cut on the tt~tt° mass is 
removed. For events in the B signal region (\Mb — 5.280 1 < 0.006 GeV) the tt~tt° spectrum is 
examined after subtracting the tt~tt° spectrum from the low B mass sideband (5.2 < Mb < 
5.26 GeV). The background subtracted tt~tt° invariant mass spectrum is then fitted to the 
sum of a Breit Wigner and a polynomial parameterization of non-resonant B° — > D* + tt~tt° 
obtained from a Monte Carlo simulation. Fig. 11 shows the fit to the background subtracted 
7T~7r° invariant mass spectrum. The fit gives an upper limit of less than 6 non-resonant B° — > 
D* + TT~Tr° events in the p mass window at the 90% confidence level. This implies that the 
non-resonant contribution to the B° — > D* + p~ decay is less than 9% at the 90% confidence 
level. If we instead take the shape of the non-resonant component from a D**(2420)7r + Monte 
Carlo we obtain a similar limit for the non-p component. A similar study has been made of 
B~ — > D*° p~ which also shows a negligible non-resonant component. The branching ratios 
for B — > D* p can be found in Tables V and VI. In Fig. 12 we show the Mb distributions 
for B~ — > D* 7r~7T~7r + and B° — > D* + TT~TT~Tr + where the 7r~7r~7r + invariant mass is required 
to be in the interval 1.0 < 7r~7r~7r + < 1.6 GeV. To show that this signal arises dominantly 
from a~[ we also present the Mb distributions for the ai sidebands 0.6 < 7r~7r~7r + < 0.9 GeV 
and 1.7 < 7r~7r~7r + < 2.0 GeV (Fig. 13), where there are signals of 15 ± 6 and ± 5.5 events 
for the D* + and D*° channels respectively. The sideband signals are 18±6% (0±13 %) of 
the signals in the ai peak, as compared to the expectation of about 10% from the tails of a 
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Breit-Wigner distribution. In Figs. 14 and 15 we show the 7r~7r~7r + mass distributions for 
aB-> D*a^ Monte Carlo simulation, a5-> D*tt~ p° non-resonant background simulation, 
and the data events in the B signal region with the scaled B mass sideband subtracted. 
The a\ meson is parameterized in the Monte Carlo simulation as a Breit-Wigner resonance 
shape with m ai = 1182 MeV and T a i = 466 MeV. The fit gives upper limits of less than 4.2 
and 4.6 non-resonant events at the 90% confidence level. This implies that the non-resonant 
components in this decay are less than 9.4% and 10.6% at the 90 % confidence level. We 
have verified that a D**(2420) / o~ Monte Carlo simulation gives a similar limit for the non-ai 
component. Our results for B meson decays into final states containing a D* meson and 
three charged pions are also listed in Tables V and VI. 



B. Polarization in B D* + p decays 



The sample of fully reconstructed B — > D* + p~ decays can be used to measure the D* + 
and p~ polarizations. By comparing the measured polarizations in B° — > D* + p~ with 
the expectation from the corresponding semileptonic B decay a test of the factorization 
hypothesis can be performed (see section XC). 

The polarization is obtained from the helicity angle distribution. The p helicity angle © p 
was defined earlier. The D* + helicity angle is the angle between the tt + direction and 
B direction in the D* + rest frame. 

The momentum in the laboratory for pions from the D* + decay which are emitted in 
the backward hemisphere (cos < in our convention) extends from 160 MeV/c down to 
about 100 MeV/c. In this momentum range, the reconstruction efficiency for charged tracks 
is reduced and becomes momentum dependent. 

Before examining the B° — > D* + p~ decay mode, we perform a consistency check of the 
efficiency correction and analysis procedure by measuring the polarization in B° — > D* + tt~. 
Since B mesons and pions are pseudoscalars, the D* + mesons from the decay B° — > D* + tt~ 
will be longitudinally polarized giving a cos 2 distribution. The same procedure used 
in the analysis of the B° — > D* + p~ polarization is applied to this case. After performing 
the sideband subtraction and correcting for efficiency, [27] we obtain the D* + helicity angle 
distribution shown in Fig. 16 (c). A fit to this distribution gives T^/T = 106 ± 7% which is 
consistent with the expectation from angular momentum conservation of T^/T = 100%. 

We now proceed to measure the polarization in B° — > D* + p~ decays. After integration 
over x, the angle between the normals to the D* + and the p~ decay planes, the helicity angle 
distribution can be expressed as follows [28]: 

— oc isin 2 D .sin 2 O p (|# +1 | 2 + |#_! | 2 ) + cos 2 © D , cos 2 © p |# | 2 (7) 

a cos fc)£>*acos fc) p 4 

The fraction of longitudinal polarization is defined by [28] 

^ L 



T r, I Hn 1 2 



r |#+i| 2 + |#-i| 2 + |#o| 2 



f8) 



If longitudinal polarization dominates, both the D* + and the p~ helicity angles will follow 
a cos 2 distribution, whereas in the case of transverse polarization we will observe a sin 2 
distribution for both helicity angles. 
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To measure the polarization we combine the helicity angle distributions for the three D° 
submodes in the B signal region (defined by \Mb — 5.280 1 < 0.006 GeV) and then subtract 
the helicity angle distribution of the scaled sideband (defined by 5.200 < Mb < 5.260 GeV). 
We fit the resulting helicity angle distributions to the functional form given in equation (6). 

From the fit to the D* + helicity angle distribution, we find T^/T = (85 ± 8)%, and from 
the corresponding fit to the p helicity angle distribution we find T^/T = (97 ±8)%. The 
results of the fits [29] are shown in Fig. 16(a) and (b) . The statistical error can be reduced 
by taking advantage of the correlation between the two helicity angles (See Fig. 17). The 
most precise result can be extracted by performing an unbinned two dimensional likelihood 
fit to the joint (cos , cos P ) distribution. This method gives 



The systematic error contains the uncertainties due to the background parameterization and 
the detector acceptance. 



In addition to the production of D and D* mesons, the charm quark and spectator anti- 
quark can also hadronize as a D** meson. The D**°(2460) has been observed experimentally 
and identified as the J p = 2 + state, while the D**°(2420) has been identified as the 1 + state 
[14]. These states have full widths of approximately 20 MeV. Two other states, a + and 
another 1 + are predicted but have not yet been observed. Presumably this is due to their 
large intrinsic widths. There is evidence for D** production in semileptonic B decays [30,31], 
and it is possible that the D** can also be seen in hadronic B decays. 

In order to search for D** mesons from B decays we first study the final states B~ — > 
D* + tt~tt~ and B~ — > D* + 7r~7r~7r°. In the latter case we require that one tt~tt° invariant 
mass is consistent with the p~ mass. The reactions B~ — > D* + tt~tt~ and B~ — > D* + 7r~7r~7r° 
have not been observed clearly in past experiments [3,6] and are not expected to occur in a 
simple picture in which the c quark plus spectator antiquark form a D* . We combine the 
D* + with a 7r~ to form a D** candidate. D** candidates lying within one full width of the 
nominal mass of either a D**°(2420) or a D**°(2460) are then combined with a tt~ or p~ to 
form a B~ candidate. 

We have also searched for D** production in the channels D + tt~tt~ and D tt~tt + . Since 
£>**°(2420) -> Dtt is forbidden, we only search for £>**°(2460) in the Dtttt final state. For 
this subset of modes, we require the Dtt mass to lie within ± 1.5 T (±28 MeV) of the nominal 
£>**(2460) mass. 

Figs. 18 and 21 show the B mass distributions for combinations of D**°(2460) or 
D**°(2420), and tt~ or p~ . In the D**°(2420)7r~ mode, there is an excess of 8.5 events 
in the B peak region with an estimated background of 1.5 events. The probability that the 
excess is due to a background fluctuation is 4 X 10~ 6 which indicates that this is a signif- 
icant signal. In this channel we give the branching ratio in Table VII, while for the other 
five combinations where the probability that the observed events are the result of a back- 
ground fluctuation is larger, we give upper limits. Our results are consistent with theoretical 
predictions [32,33] based on the factorization hypothesis (Table VIII). 




(9) 



VI. MEASUREMENTS OF D** FINAL STATES 
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We have also investigated the final states D* + tt~tt~ (D + tt~tt~ and D tt + tt~) where the 
and a charged pion are not constrained to lie in any particular mass interval. To observe 
these signals, the background from the final states p~ must be suppressed. Since the 
p~ final state is highly polarized (see section VB), it is possible that the slow tt° from the 
p~ decay can be exchanged for a slow charged pion from the decay of the other B meson. To 
eliminate this background, we make cuts on the cosines of the helicity angles, and ©£>**, 
where &£>** is calculated for the D* + tt~ 1ow system with 7r s ^ ow being the slower tt~ of the two. 
This helicity angle is defined as the angle between the B and the fast tt~ in the rest frame of 
D* + 7r s ^ ow system. We require cos &£>** < 0.8 and | cos ®r>* \ < 0.7. For the B — > Dtttt modes, 
a similar cut, cos P < 0.7, is made using the tttt system. In Figs. 20 and 19, we show the 
Mb candidate mass distributions. There is a significant signal in B — > D* + tt~tt~ . For the 
other two modes we quote upper limits in Table VII. 

VII. EXCLUSIVE B -> CHARMONIUM DECAYS 
A. Introduction 

In B decays to charmonium the c quark from the b combines with a c quark from the 
virtual W~ to form a charmonium state. This process is described by the color suppressed 
diagram shown in Fig. 1(b). By comparing B meson decays to different final states with 
charmonium mesons the dynamics of this decay mechanism can be investigated. The decay 
modes B° — > ipK° and B° — > ip'K are of special interest since the final states are CP 
eigenstates which can be used to determine one of the three CP violating angles accessible 
to study in B decays. It is also possible to use the B° — > ij)K*° decay (where K*° — > Kgir ) 
to measure this CP asymmetry. However, this final state has even CP if the orbital angular 
momentum L, between the ip and K*° is or 2, and odd CP for L = 1. If both CP states are 
present the CP asymmetry will be diluted. A measurement of CP violation in this channel 
may be possible if one of the CP states dominates, or if a detailed moments analysis of the 
various decay components is performed [34]. We present a measurement of the polarization 
in the decay B° — > if)K*° which allows us to determine the fractions of the two CP states. 

B. Branching Ratios 

B meson candidates are formed by combining a charmonium and a strange meson can- 
didate. We reconstruct Kg decays into tt + tt~ pairs which have vertices displaced from the 
beam axis by greater than 5 mm. Only the K + tt~ channel is used to form K*° mesons while 
K*~ candidates are reconstructed in the decay channels KgTV~ and K~tt°. The Ktt invariant 
mass must be within ±75 MeV of the nominal K* mass. In the B~ — > ipK*~ , K*~ — > K~tt° 
mode, only the half of the K*~ helicity angle distribution with a fast 7r° is used. The ij) and 
%\)' candidates decaying into lepton pairs were kinematically constrained to the known mass 
values in order to improve the resolution on AE. Using the procedures described in section 
III we reconstruct B meson candidates and obtain the beam constrained mass distributions 
shown in Figs. 22, 23 and 24. The corresponding branching ratios are listed in Table IX. 
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The systematic errors in the branching ratio measurements include contributions from 
number of B mesons (2.5%), tracking efficiencies (2% per charged track), 7r° detection effi- 
ciency (5%), dE I dx efficiency (2% per identified track), lepton detection efficiency (2% per 
lepton), Monte Carlo statistics (1.5 - 6%), the ip leptonic branching ratio (4.2%) and the 
branching ratios for ip and Xci decays. 

The results for the B~ and B° decay modes can be combined using isospin symmetry to 
determine the vector to pseudoscalar production ratio 



Bill rl\ 
Bill -rlxi 



1.71 (10) 



The revised Bauer-Stech-Wirbel (BSW) model [8] predicts a value of 1.61 for this quantity. 
This model uses the ratio of B — > K* / B — > K form factors determined from harmonic 
oscillator wavefunctions and assumes that the factorization hypothesis is valid for internal 
spectator decays. 



C. Polarization in ipK* 

After integration over the azimuthal angle between the ip and the K* decay planes, the 
angular distribution in B — > ipK* decays can be written as [28] 

d 2 T 1 

- — — — oc - sin 2 0^.(1 + cos 2 0^)(|# +1 | 2 + |#_i| 2 ) + cos 2 A -» sin 2 0^|# o | 2 , 

d cos Ky^pd cos k)k* 4 

(11) 

where the K* helicity angle ®k* is the angle between the kaon direction in the K* rest frame 
and the K* direction in the B rest frame, 0^, is the corresponding ip helicity angle, and the 
H±i j0 are the helicity amplitudes. 

There are 29 B° — > if)K*° candidates and 13 B~ — > ipK*~ candidates. After correcting 
for detector acceptance, we perform an unbinned maximum likelihood fit to the double differ- 
ential distribution described in the equation above. The fit gives the fraction of longitudinal 
polarization in B — > ipK* as 

(^f) =0.80 ±0.08 ±0.05 (12) 

V 1 J B^tpK* 

The systematic error in this measurement is dominated by the uncertainty in the acceptance. 
The efficiency corrected distributions for each of the helicity angles cos 0^, and cos ®k* are 
shown in Fig. 25. 

This result can be compared to the theoretical predictions of Kramer and Palmer [66] 
which depends on the unmeasured B — > K* form factor. Using the BSW model to estimate 
the form factor, they find T^/T = 0.57. Using Heavy Quark Effective Theory (HQET) and 
experimental measurements of the D — > K* form factor, they obtain T^/T = 0.73. 

The decay mode B — > ipK* may not be completely polarized, but it is dominated by a 
single CP eigenstate (CP = —1 produced with L = 1). This mode will therefore be useful 
for measurements of CP violation. 
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VIII. SEARCH FOR COLOR SUPPRESSED DECAYS 



We search for B decays which can occur via an internal VK-emission graph, but which 
do not lead to final states with charmonium [39]. One expects that these decays will be 
suppressed relative to decays which occur via the external VK-emission graph. For the internal 
graph the colors of the quarks from the virtual W must match the colors of the c quark and 
the accompanying spectator antiquark. In a simple picture, one expects that the suppression 
factor should be about 1/18 for decays involving 7r°, p° and u mesons [36], but in heavy quark 
decays the effects of gluons cannot be neglected. These decays can be used to test QCD 
based calculations [8] which predict suppression factors of order 1/50. If color-suppressed B 
decay modes are not greatly suppressed then these modes could be useful for CP violation 
studies [38]. 

We search for color-suppressed decay modes of B mesons which contain a single D meson 
(or D* meson) in the final state. The relevant color-suppressed modes are given in Table X. 
We use the decay modes rj — > ~{~{ } to — > 7r + 7r~7r and rj — > ?y7r + 7r~, followed by ?y — s- 77 [37]. 

For decays of a pseudoscalar meson into a final state containing a pseudoscalar and a 
vector meson, a helicity angle cut of |cos0y| > 0.4 is used [40]. No convincing signals 
were found in the decay modes that were examined. Upper limits on the branching ratios 
for color-suppressed modes are given in Table X. The 90% confidence level upper limits 
are calculated using the prescription described by the PDG [41]. These upper limits take 
into account the systematic uncertainty in the background level as well as the systematic 
uncertainty in the detection efficiency. In Figs. 26, 27, and 28 we show the fitted distributions 
for the color-suppressed modes with the fit superimposed on each plot. Upper limits on the 
ratios of color-suppressed modes to normalization modes are given in Table XL 

IX. THE B~ - B° MASS DIFFERENCE 

We now proceed to measurements of the B° and B~ masses and the mass difference 
between them. For this analysis we use the decays B~ — > ipK~, B° — > ipK*°, B~ — > D°tt~ } 
B~ -> D°p~, B~ -> D*°tt~, B~ -> D*°p~, B° -> D+tt", B° -> D+ p~ , B° -> D* + Tr~, and 
B° — > D* + p~ for which the signal to background ratio is large. For the decays B~ — > D°tt~ 
and B~ — > D°p~ only the D° — > K~tt + mode is used. The Mb distributions for the sum 
of these modes are shown in Fig. 29. We have a total of 362 B~ and 340 B° signal events. 
The data are fitted with a Gaussian of fixed width (2.7 MeV) determined by Monte Carlo 
simulation. The width is assumed to be the same for all modes. The fitted masses for each 
mode and their statistical errors are given in Table XII. We apply a correction for initial 
state radiation as described in Ref. [42], of magnitude - 1 . 1 ±0.5 MeV, to arrive at the final 
values for the B~ and B° masses of 5278.8 ± 0.2 ± 0.5 ± 2.0 MeV and 5279.2 ± 0.2 ± 0.5 ± 2.0 
MeV, respectively. The first systematic error results from the uncertainty in the initial state 
radiation correction. The second systematic error is due to the uncertainty in the absolute 
value of the CESR energy scale, which is determined by calibrating to the known T(1S) 
mass [43]. 

The mass difference is determined to be 0.41±0.25±0.19 MeV. This is more accurate than 
the masses themselves because the beam energy uncertainty cancels, as do many systematic 
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errors associated with the measurement errors on the charged tracks and 7r° mesons. The 
remaining systematic error is found by making a number of tests of the stability of our result. 

A systematic shift of 0.12 MeV is produced by using different background shapes for the 
B~ and B° modes [44]. We have also investigated the effect of changing the photon energy 
calibration. A change of 0.5%, the quoted systematic error, results in a 0.15 MeV change in 
the fitted B mass in both the D* + p~ and D*°p~ final states. This effect almost completely 
cancels in the mass difference measurement where it contributes an error of <0.03 MeV. 
This is because the shift of the 7r° energy in the p~ cancels in the difference leaving only the 
shift of the energy of the slow 7r° from the D*° which is uncorrelated with the direction of 
the B meson. A similar test where we scale the measured momentum of the slow pion from 
the D* + by 5 %, also does not affect the mass difference for the same reason. 

We have also checked the stability of the result with changes in event samples. For 
example, we have used only half of the cos P distributions in B — > p modes and we 
have used less stringent lepton identification criteria for the B — > ipK~ mode. We estimate 
a systematic of 0.15 MeV from these studies. 

The different sources of systematic errors are listed in Table XIII and are combined in 
quadrature. We compare our result with previous results in Table XIV. 

There are several models which predict the isospin mass difference to be between 1.2 
and 2.3 MeV which are larger than the value reported here [45]. However, Goity and Hou 
(—0.5 ± 0.6 MeV) and Lebed (0.89 MeV) can accomodate this small mass difference in their 
models [46]. 

X. TESTS OF THE FACTORIZATION HYPOTHESIS 

A. Introduction 

Our large data sample has made possible the precise branching ratio and polarization 
measurements discussed above. In the following sections we address many important ques- 
tions about non-leptonic B meson decay. 

By comparing rates and polarizations of semileptonic and hadronic decays we can perform 
tests of the factorization hypothesis, which is the basis of most theoretical treatment of 
hadronic B decays. In analogy to semileptonic decays, where the amplitude factorizes into 
the product of a leptonic and hadronic current since leptons are not sensitive to the strong 
interaction, it is possible that two body decays of B mesons which occur via the external 
spectator process may be expressed theoretically as the product of two independent hadronic 
currents, one describing the formation of a charm meson and the other the hadronization of 
the ud (or cs) system from the virtual W~ . 

There are few models of hadronic B decays. Those which exist predict widths of two- 
body decays and assume the validity of the factorization hypothesis. Although factorization 
fails in many D decays [47], it is hoped that factorization will be a better approximation in 
B decays due to the larger energy release present [48]. 

If factorization is valid, then heavy quark effective theory, henceforth referred to as HQET 
[49], could provide a reliable, model independent framework for the calculation of properties 
of non-leptonic B meson decays. In addition, if factorization holds, then measurements 
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of non-leptonic B decays may be used to extract fundamental parameters of the Standard 
Model. For instance the CKM matrix element V ub can be determined from B° — > tt + tt~ or 
B° — > D~tt + , and the decay constant fr> s can be determined from B° — > D~ D* + . 



B. Branching Ratio Tests 

Assuming factorization, the effective Hamiltonian Eq. (2) for a non-leptonic B decay can 
be written as a product of two hadronic currents. Consider the case of B° — > D* + h~ , where 
h is a hadron. The amplitude for this reaction is 

A = G F /V2 V cb V: d (h-(p)\(du)\0)(D* + \(cb)\B^) (13) 

where V u d is the well measured CKM factor from the W~ — > ud vertex. The first hadron 

current, which creates the h~ from the vacuum, is related to the decay constant and is 
known for h = 7r, p. We have 

(h-(p)\(du)\0) = -if h p,, (14) 

where p^ is the h~ four momentum. The other hadron current can be found from semileptonic 
B° — > D* + £~ ' vi decays. Here the amplitude is the product of a lepton current and the hadron 
current that we seek to insert in Eq. (13). Factorization can be tested experimentally by 
verifying whether the relation 



r B° -> D*+h- 



rIV 



6x'4n\V^\', (15) 



is satisfied, where q 2 is the four momentum transfer from the B meson to the D* meson. 
Since q 2 is also the mass of the lepton-neutrino system, by setting q 2 = we are simply 
requiring that the lepton-neutrino system has the same kinematic properties as the h~ in the 
hadronic decay. The c\ term accounts for hard gluon corrections. Here we use Ci = 1.1 ± 0.1 
as deduced from perturbative QCD. The error in Ci reflects the uncertainty in the mass scale 
at which the coefficent Ci should be evaluated [50]. For the case where h~ = tt~ and Ci = l, 
equation 15 was found to be satisfied by Bortoletto and Stone [51]. In the following the left 
hand side of Eq. (15) will be denoted i? exp and the right hand side will be denoted i?Theo- 

This type of factorization test can also be performed using B° — > D* + h~ decays where 
h~ = p~ or ai(1260)~. For the p~ case Eq. (15) becomes: 

r = r(ff°^D*V) = 2 2 2| |2 ( } 

§ (B ^ D*l v)\ qa=m% lJpl 1 ; 

where the semileptonic decay is evaluated at q 2 = m 2 = 0.60 GeV 2 . The decay constant on 
the right hand side of this equation can be determined from e + e~ — > p° or from r decays. 
The first method leads to f p = 215 ± 4 MeV. Taking into account the p width, Pham and 
Vu [52] find that T(t~ -> vp~) = 0.804(^/16^ \V 2 d \M^f 2 which gives f p = 212.0 ± 5.3 
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MeV [53]. We take the first value. We also perform this test for B° — > D* + a~[ where we 
use f ai = 205 ± 16 MeV [33]. To derive numerical predictions for branching ratios, we must 
interpolate the observed differential q 2 distribution [54] for B — > D*£ v to q 2 = m 2 } m 2 } 
and m 2 , respectively. Until this distribution is measured more precisely theoretical models 
must be used for the slope of the distribution. Thus the results are stated below for different 
models. Fortunately, the spread in the theoretical models which describe B — > D*£ v is small 
(see Fig. 30). 

We now have all the required ingredients [55] for the test with decay rates (see Table 
XV). Using the extrapolation of the q 2 spectrum [56] from the WSB model as the central 
value, we obtain from Eqs. (15) and (16) the results given in Table XVI. 

If we form ratios of branching fractions some of the systematic uncertainties on R exp will 
cancel, as does the QCD correction c\ in Rtheor- For example in the case of D* + p~ / D* + tt~ , the 
expectation from factorization is given by Rtheor (/?)/ Rtheor (tt) times the ratio of the semilep- 
tonic branching ratios evaluated at the appropriate q 2 values. In Table XVII we show the 
comparison of the data, the expectation from factorization as defined above and two the- 
oretical predictions of Bauer, Stech and Wirbel (BSW) [64], and Reader and Isgur (RI) 
[33]. 

From the measurements described above, we find that at the present level of precision, there 
is agreement between experiment and the expectation from factorization for the q 2 range: 
< q 2 < m 2 ai . 

C. Factorization and Angular Correlations 

More subtle tests of the factorization hypothesis can be performed by examing the po- 
larization in B meson decays into two vector mesons. This idea was suggested by Korner 
and Goldstein [57]. Again, the underlying principle is that hadronic decays are analogous to 
the appropriate semileptonic decays evaluated at a fixed value of q 2 . For instance, the ratio 
of longitudinal to transverse polarization (Tl/Tt) in B° — > D* + p~ should be equal to the 
corresponding ratio for B — > D*l~ v evaluated at q 2 = m 2 = 0.6 GeV 2 . 




The advantage of this method is that it is not affected by QCD corrections [58]. 

For B — > D* I v decay, longitudinal polarization dominates at low q 2 . Near q 2 = q^ ax7 
by contrast, transverse polarization dominates. There is a simple physical argument for the 
behaviour of the form factors near these two kinematic limits. Near q 2 = ^ ax , the D* is 
almost at rest. Its small velocity is uncorrelated with the D* spin, so all three possible D* 
helicities are equally likely. As q 2 — > ^ ax we expect Yt/Yl = 2. At q 2 = 0, the D* has the 
maximum possible momentum, while the lepton and neutrino are collinear and travel in the 
direction opposite to the D* with their helicities aligned to give S z = 0. Thus, near q 2 = 
longitudinal polarization is dominant. 

For B° — > D* + p~ , Rosner predicts 88% longitudinal polarization from the argument 
described above [59]. Similar results can be extracted from the work of Neubert [60] and 
Kramer et al. [68]. Fig. 31 shows Neubert's result for the production of transversely and 
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longitudinally polarized D* mesons in B — > D* I v decays. Using this figure we find T^/T to 
be approximately 85% for q 2 = m 2 = 0.6, which agrees well with Rosner's prediction [59]. 
The agreement between these predictions and the experimental result (see section VB) 

r L /r = (93±5±5)% (18) 

supports the factorization hypothesis in hadronic B meson decay for q 2 values up to m 2 . 



XI. TESTS OF SPIN SYMMETRY IN HQET 

If the factorization hypothesis holds, then certain hadronic B meson decay modes can 
be used to test the spin symmetry of HQET. In this theory the effect of the heavy quark 
magnetic moment does not enter to lowest order [61] , so it is expected that 

T(B° -> D + Tr~) = T(B° -> D* + ir~) (19) 



anc 



T(B° -> D+p~) = T(B° -> D* + p~). (20) 

After correcting for phase space and deviations from heavy quark symmetry, one expects 
that B(B° -> D+TT-) = 1.03 B(B° -> D*+ir-) and B(B° -> D+p~) = 0.89 B(B° -> D* + p~). 
A separate calculation by Blok and Shifman using a QCD sum rule approach predicts that 
B(B° — > D + tt~) = 1.2 B(B° — > D* + tt~) due to the presence of non-factorizable contributions 
[62]. From our data we find 



B(B° -> D+ 



7T 



B{B° -»■ D*+tt- 

and 

B(B° -> D+p- 



B(B° -> D*+p~ 



1.12 ±0.19 ±0.24 (211 



1.10 ± 0.14 ± 0.28 (22) 



The contribution in this ratio from the systematic error on the detection efficiency is re- 
duced to 5% for these two cases. Both ratios of branching fractions are consistent with the 
expectation from HQET spin symmetry as well as the prediction from Blok and Shifman 
[62]. 

Mannel et al. [61], also observe that by using a combination of HQET, factorization, and 
data on B — > D*£v from CLEO and ARGUS they can obtain model dependent predictions 
for B(B° — > D + p~) I B(B° — > D + tt~). With three different parameterizations of the B — >■ D 
form factor [63] this ratio is predicted to be 3.05, 2.52, or 2.61. 

From the measurements of the branching ratios we obtain 

B(B°^D + p-) , , 

— H ^— f = 2.8 ±0.5 ±0.2 23 

B{B° -> D+TT-) y J 

The systematic errors from the D branching fractions and the tracking efficiency cancel in 
this ratio. Thus we find good agreement with the prediction from HQET combined with 
factorization. 



19 



XII. DIRECT AND INDIRECT EFFECTS OF THE COLOR SUPPRESSED 

AMPLITUDE 



A. Introduction 



In the QCD treatment described by equations (1) and (2) it is difficult to take into account 
the effects of multiple soft gluon emission analytically. Instead, in the phenomenological BSW 
approach [64] two undetermined coefficients are assigned to the effective charged current, 
ai(/i), and the effective neutral current, a 2 (fi) } parts of the B decay Hamiltonian. These 
coefficients were determined from a fit to a subset of the experimental data on charm decays. 
With these values the decay rates for a large number of non-leptonic decays can then be 
calculated using the factorization hypothesis, and model dependent hadron form factors. 
We can relate ai(/i) and a 2 (fi) to the QCD coefficients Ci(/i) and c 2 (/i) by ai = Ci + ^c 2 
and a 2 = c 2 + <^Ci where £ = l/iV co i or . The values ai(m 2 c ) = 1.3 and a 2 (m 2 c ) = —0.55 which 
give the best fit to the experimental data on charm decay correspond to l/iV co i or ~ [8]. 
However, there is no rigorous theoretical justification for this choice of iV co i or [65]. 

In the decays of charmed mesons the effect of color suppression is obscured by the effects of 
final state interactions (FSI) and soft gluon effects which enhance W exchange diagrams. For 
instance, Table XVIII gives ratios of several charmed meson decay modes with approximately 
equal phase space factors where the mode in the numerator is color suppressed while the 
mode in the denominator is an external spectator decay [67]. These modes are clearly not 
suppressed. However, the following decay appears to be suppressed. 

B(D° -> K°p°) 

[ P ' 0.08 ± 0.04 (24) 



B(D° -> K-p+) 



In contrast to the charm sector where the mechanism of color suppression is obscured, 
one expects to find in B meson decays a simple and consistent pattern of color suppression. 
Partly, it is expected that color suppression is more effective at the b quark mass scale than 
the charm quark mass scale due to the evolution of the strong coupling constant a s to smaller 
values. Using the BSW model and extrapolating from q 2 = m 2 c to q 2 = ml using the values 
from charm decays, gives the predictions ai(ra^) = 1.1 and a 2 (ml) = —0.24 for B decays. 
Another approach using the factorization hypothesis, HQET and model dependent form 
factors has been suggested by C. Reader and N. Isgur (RI model) [33]. In this approach, a\ 
and a 2 are determined from QCD (with l/iV co i or = 1/3) and color suppressed B decays are 
expected to occur at about 1/1000 the rate of unsuppressed decays. Observation of these 
decays at a much greater level would indicate the breakdown of the factorization hypothesis. 
In section VIII we obtained upper limits for color suppressed B decays with a D° or D*° 
meson in the final state. In Table XIX these results are compared to prediction of the BSW 
and the RI model. 

In contrast to charm decays, color suppression seems to be operative in hadronic decays 
of B mesons. The limits on the color suppressed modes with _D°(*) and neutral mesons are 
still above the level expected in the model of Bauer, Stech and Wirbel. However, the limit 
on B° -> D°tt° disagrees with Terasaki's prediction [70] that B(B° -> D°tt°) w 1.8 B(B° -> 
D + tt~). To date, the only color suppressed B meson decay modes which have been observed 
are final states which contain charmonium mesons e.g. B — > ipK and B — > ipK* [71]. 
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B. Determination of |ai|, | «2 1 an d the relative sign of (a 2 /ai) 

In the BSW model [8,64] , the branching tractions of the B° normalization modes are 
proportional to a\ while the branching fractions of the B — > if) decay modes depend on 
a\ (Table XX [8]). A fit to the branching ratios that we have measured for the modes 
B° -> D+7T", D+p~, D*+tt- and D*+ p~ yields 

H = 1.15 ±0.04 ±0.05 ±0.09 (25) 

and a fit to the modes with ip mesons in the final state gives 

|a 2 | = 0.26 ±0.01 ±0.01 ± 0.02 (26) 

The first systematic error on | «i | and \a 2 \ includes the experimental uncertainties from the 
charm or charmonium branching ratios, tracking efficiency, background shapes and the value 
of \V c b\ } but does not include the theoretical uncertainties. There is a second uncertainty due 
to the B meson production fractions and lifetimes. These are constrained by the value of 
(f+T + /f T ) determined from the CLEO II [72] measurement of B(B~ -> D*° 7" v) / B (B° -> 
D* + l~v) = 1.20 ± 0.20 ± 0.19. 

The comparison of B~ and B° modes can be used to distinguish between the two possible 
choices for the sign of a 2 relative to a\. The BSW model, Ref. [8] predicts the following ratios: 

= B(B~ ^D°7r-) h 23a 2 /a 1 ) 2 (27) 

B B°^D+tt-) y 1 ' y ' 



B(B-^D°p-) , , x9 

R* = g^o^^-j = C 1 + O-eea./a,) 2 (28) 

The numerical factor which multiplies a 2 /ai is proportional to the ratio of B — > to 
B — > Tr(p) form factors as well as the ratio of the 7r(p) meson to D meson decay constants. 
We assume fr> = fp* = 220 MeV [69]. Only the B — > D* form factor and the 7r(p) meson 
decay constant have been measured experimentally. 
Similarly, we define 



B(B~ — >■ D*°tt~) , , x9 

^3 = ^ = 1 + 1.29a 2 / fll 2 (29) 



12(f)- r)*°n~\ 

o^d £L_£_i«(i + .75a 2 /ai) 2 (30) 

B(B° D*+p~) y J 

Table XXI shows a comparison between the experimental results and the two allowed solu- 
tions in the BSW model. In these ratios, the systematic errors due to detection efficiency are 
reduced. In the ratios R 3 and i? 4 the D° — > K~tt + branching ratio error does not contribute 
to the systematic error. 

It is important to note that the determination of the sign of a 2 /ai depends on assumptions 
about the relative production of B + and B° mesons at the T(4S') resonance, / + and / , as 
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well as their lifetimes, r + and r . A least squares fit to the above ratios using the CLEO II 
value for (f + T + /f T ) [72] gives a 2 /ai = 0.23 ± 0.04 ± 0.04 ± 0.10 where we have ignored 
uncertainties in the theoretical predictions for Ri through i? 4 . The second systematic error 
is due to the uncertainty in (f + T + /f T ). As this ratio increases, the value of a 2 /ai decreases. 
The allowed range of (f + T + /f T ) excludes a negative value of a 2 /ai. Other uncertainties in 
the magnitude of fr> and the B — > tt form factor can change the magnitude of a 2 /ai but not 
its sign. This result is consistent with the value of a 2 determined from the fit to the B — >■ if) 
decay modes. It disagrees with the theoretical extrapolation from data on charmed meson 
decay in the BSW model [73] which predicts a negative value for a 2 /ai. 

XIII. CONCLUSIONS 

We have presented new measurements of B branching ratios, resonant substructure and 
masses. More accurate branching ratios are given for many modes. The modes B — > D* p~ 
and B — > D*a~[ are clearly seen for the first time. 

Using a subset of 702 B meson decays reconstructed in channels with good signal to 
background ratios we have made a precise measurement of the B° — B~ mass difference of 
0.41 ± 0.25 ± 0.19 MeV. 

We have carried out an extensive series of tests of the factorization hypothesis including 
comparisons of rates for D* + h~ (where h~ = tt~ } p~ } or a~[) with rates for D* + l~v at q 2 = 
M|, as well as comparisons of the polarizations in D* + p~ with D* + £~Pi. In all cases the 
factorization hypothesis is consistent with the data. 

We have made improved measurements of branching ratios of two-body decays with a 
if) } if) or Xc meson in the final state. The decay B — > if)K* is strongly polarized with 
Tl/T = 0.80 ± 0.06 ± 0.08. Therefore this mode will be useful for measuring CP violation. 

A search for color suppressed decays with a charmed meson and light neutral hadron 
in the final state shows no positive evidence for such processes. The most stringent limit, 
B(B° — > D°tt°) I B(B° — > D + tt~) < 0.09, is still above the level where these color suppressed 
B decays are expected in most models. 

The observation of B — > if) modes shows that color suppressed decays are present. Using 
only exclusive B — > if) decays we find a value of the BSW parameter \a 2 \ = 0.26±0.0I±0.0I± 
0.02. We also report a new value for the BSW parameter | «i | = 1. 15 ± 0.04 ± 0.05 ± 0.09. 
Comparing B + and B° decays, we find a 2 /ai = 0.23 ± 0.04 ± 0.04 ± 0.10. We have shown 
that the sign of a 2 /ai is positive, in contrast to what is found in charm decays. 
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TABLES 



TABLE I. D° and D+ decay modes 

D type Decay Mode B(%) a mo (MeV) 

D 6 K~tt+ 3.91±0.08±0.17 8L5 

D° K+tt-tt 12.1±1.1 13.0 

D° K-tt+tt+tt- 8.0±0.5 8.1 

D+ K-ir + ir + 9.1 ±1.4 7.6 



TABLE II. D* decay modes used 

D* Mode B{%) cr mct _ mc (MeV) 

D*+ D°ir+ 68.1±1.6 0.8 

D*° D°ir° 63.6±4.0 1.1 



TABLE III. Branching Ratios (%) for B~ D°(mr)- 



B~ Mode 


D mode 


& AE 


#of 




B(%) 


B average (%) 






(MeV) 


events 










K~ir + 


22 


76.3 ± 9.1 


0.433 


0.48±0.06 






K-ir+ir 


26 


134 ± 15 


0.193 


0.62±0.07 


0.55±0.04±0.05± 0.02 




K~ir + ir + 'ir~ 


20 


94± 11 


0.222 


0.57±0.07 






K~ir + 


18-38 


80 ± 9 


0.155 


1.40 ± 0.18 




D°p- 


K-ir+ir° 


22-42 


42 ± 9 


0.036 


1.04±0.23 


1.35±0.12±0.14± 0.04 






17-37 


90.4±12.1 


0.079 


1.53±0.20 




a This efficiency does not include D branching ratios. 








TABLE IV. Branching Ratios 


(%) for B° -+ D+(mr) 


r 


B° Mode 


D Mode 


& AE 


#of 




B(%) 


B average (%) 






(MeV) 


events 








D+ir- 




20.5 


80.6±9.8 


0.32 


0.29±0.04 


0.29±0.04±0.03± 0.05 


D+p- 




18-38 


78.9±10.7 


0.12 


0.81±0.11 


0.81±0.11±0.12± 0.13 



This efficiency does not include D branching ratios. 
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TABLE V. Branching Ratios (%) for B~ D*°(mr) 



B~ Mode 


D° 


Mode 


& AE 

(MeV) 


#of 

events 


e a 


B(%) 


B average (%) 




K~ 


■7T + 


25 


13.3±3.8 


0.16 


0.36±0.13 




D 7T 


A 


IT IT 


oz 


o < . i ±o.y 


U.Uo 


n coin 1 o 


U.OZ±U.U f ±U.Ub ± U.U4 




K~ 


1T + 1T + 1T~ 


21 


20.0±4.9 


0.08 


0.52±0.13 






K~ 


TT + 


21-41 


25.7±5.4 


0.064 


1.74±0.37 




D*°p~ 


K~ 


-t- n 


26-46 


43.8±7.8 


0.027 


2.24±0.40 


1.68±0.21±0.25 ± 0.12 




K~ 


1T + 1T + 1T~ 


19-39 


16.9±4.6 


0.030 


1.19±0.35 






K~ 


TT + 


14 


5.5±2.9 


0.048 


0.51±0.26 




J D*°7r-7r-7r +b 


K~ 


1T + 1T 


22 


27.7±7.2 


0.022 


1.74±0.45 


0.94±0.20± 0.16 ± 0.06 




K~ 


'7r + 7r + 7r~ 


15 


15.0±4.5 


0.025 


1.26 ±0.37 





a This efficiency does not include D or D* branching ratios. 

b The three pion mass is required to be between 1.0 GeV and 1.6 GeV consistent with an a\ meson. 
(If this channel is dominated by aT/, the branching ratio for D*°a^ is twice that for _D*°7r _ 7r~7r + .) 



TABLE VI. Branching Ratios (%) for B° D* + {mr)- 



B° Mode 


D° 


Mode 


<J AE 

(MeV) 


#of 

events 




B(%) 


B average (%) 




K~ 


■7T + 


25 


19.4±4.5 


0.35 


0.22±0.05 




D*+ir- 


K~ 


■7T + 7T 


32 


31.9±6.4 


0.14 


0.30±0.06 


0.26±0.03±0.04± 0.01 




K~ 


'7r + 7r + 7r~ 


21 


20.5±5.2. 


0.15 


0.27±0.07 






K~ 


■7T + 


21.5-41.5 


21.9±5.2 


0.12 


0.71±0.17 




D*+p~ 


K~ 


■7T + 7T 


23-43 


39.8±7.2 


0.048 


1.08±0.20 


0.74±0.10±0.14± 0.03 




K~ 


'7r + 7r + 7r~ 


20.5-40.5 


14.6±4.6 


0.054 


0.52±0.17 






K~ 


■7T + 


14 


13.5±3.9 


0.096 


0.58±0.17 




D*+ir-ir-ir +h 


K~ 


■7T + 7T 


22 


21.7±5.9 


0.043 


0.67±0.18 


0.63±0.10±0.11± 0.02 




K~ 


'7r + 7r + 7r~ 


15 


13.9±4.4 


0.042 


0.65±0.19 





a This efficiency does not include D or D* branching ratios. 

b The three pion mass is required to be between 1.0 GeV and 1.6 GeV consistent with an a\ meson. 
(If this channel dominated by aT/, the branching ratio for D* + a^ is twice that for _D* + 7r~7r~7r + .) 
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TABLE VII. Branching Ratios (%) for B D**(nTr) 



B Mode 


D Mode 


CAE 




# of £> average (%) 






(MeV) 




events 


U TT TT 


Ii TT+~ 


1 7 
1 1 


n 1 q 


< 1U.1 < U.1D 


D + TT~TT~ 


K~TT + TT + 


15.5 


0.11 


< 10.3 < 0.14 


J D**(2460)vr- D+tt'tt- 


K~TT + TT + 


16 


0.21 


< 5.6 < 0.13 


j D**(2460)tt- D®TT + TT~ 


K~TT + 


17 


0.26 


< 5.6 < 0.22 


A* (24t)0J/O — ► TT TT TT 


A TT^TT^ 


lb 


0.08 


< 6.1 < 0.4/ 


D**{2Am)p- -> D°TT + TT-TT° 


K~TT + 


17 


0.11 


< 5.1 < 0.49 




K~TT + 


16 


0.161 




U TT TT 


A TT TT 


16 


U.Ubl 


1/1 1 _l_ £/l n in l n l n no l n m 
14.1 ± 0.4 U.l» ± U.lw ± U.Uo ± U.U1 




K~TT + TT~TT + 


19 


0.075 






K~TT + 


16 


0.161 




A* ^2420J7T — ► A T 7T TT 


A TT TT 


23 


0.061 


8.5 ± 3.8 0.11 ± 0.05 ± 0.02 ± 0.01 




K~TT + TT~TT + 


19 


0.075 






K~TT + 


16 


0.161 




A (24o0j7T —> JJ T 7T 7T 


A TT TT 


26 


0.061 


3.5 ± 2.3 < 0.28 




K~TT + TT~TT + 


19 


0.075 






K~TT + 


30 


0.078 




Z>**(2420)/o _ D^TT-TT-Tr 


K-TT + TT 


24 


0.037 


3.4 ± 2.1 < 0.14 




K~TT + TT~TT + 


27 


0.042 






K~TT + 


30 


0.078 




A>**(2460>- D^TT-TT-Tr 


K-TT + TT 


24 


0.037 


3.2 ± 2.4 < 0.5 




A"~7r + 7r + 7r~ 


27 


0.042 




a The efficiencies do not include the branching 


I ratios for D, D* and D**. To determine the B decay 


branching ratios, we assumed £>(_D**°(2420) - 


D* + TT- 


") and £(a>**°(2460) D* + tt~) are 67% and 


20% respectively. We also assume that B(D 


**°(2460; 


)^D+ 


tt~) and £(a>**+(2460) D°tt + ) are 



30% and 30% respectively. 
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TABLE VIII. Branching Ratio for B D* 



Mode CLEO II Bari model [32] RI model [33] 

J D**°(2420)vr- (11 ± 5 ± 2 ± 1) x 10" 4 4 x 10" 4 7.5 X 10" 4 - 13 X 10" 4 

J D**°(2460)vr- < 2.8 X 10" 3 6 X 10" 4 5 X 10" 4 - 8 X 10" 4 
(D**° D*+tt-) 

J D**°(2460)vr- < 1.3 X 10" 3 6 X 10" 4 5 X 10" 4 - 8 X 10" 4 
(D**° -> D+tt-) 

J D**+(2460)vr- < 2.2 x 10" 3 6 X 10" 4 5 X 10" 4 - 8 X 10" 4 
(D**+ -> D°ir+) 

D**°(2420)p- < 1.4 X 10" 3 1 X 10" 3 13 X 10" 4 - 24 x 10" 4 

J D**°(2460>" < 5 x 10" 3 1 x 10" 3 10 x 10" 4 - 20 x 10" 4 



TABLE IX. Exclusive B — ► cc Branching Ratios and 90% Confidence Level Upper Limits (%). 



B Mode 


a(AE) 


# of events 


e a 


B(%) 


B~ ipK~ 


13 


58.7 ± 7.9 


0.47 


0.110 ± 0.015 ± 0.009 


B° tjjK 


13 


10.0 ±3.2 


0.34 


0.075 ± 0.024 ± 0.008 


B° ipK*° 


12 


29.0 ±5.4 


0.23 


0.169 ± 0.031 ± 0.018 


H -> ipK*~, K*~ -> A -" 


21 


6.0 ± 2.4 


0.07 


0.218 ± 0.089 ± 0.026 


B~ -> ipK*~, K*~ -> K° s ir- 


11 


6.6 ± 2.7 


0.17 


0.130 ± 0.058 ± 0.018 


B~ — ► ipK*~ (combined) 




12.6 ±3.6 




0.178 ± 0.051 ± 0.023 



B~ i>'K~ 


9.8,11 


7.0 ± 2.6 


0.36,0.15 


0.061 ± 0.023 ± 0.009 


B° 4,'K 


8.4,10 





0.28,0.11 


< 0.08 




9.7,10 


4.2 ± 2.3 


0.24,0.091 


< 0.19 


H -> ip'K*~, K*~ -> A -" 


18,17 


1 ± 1 


0.077,0.023 


< 0.56 


B~ -> ip'K*~, K*~ -> K° s ir- 


7.9,9.8 


1 ± 1 


0.16,0.057 


< 0.36 


B~ — ► ip K*~ (combined) 




2 ± 1.4 




< 0.30 



B- 



B- 



B~ -> XdK- 
B° - XdK° 
B° -+ XciK* 

\ |A • A - 



K-ir° 
K s ir~ 



XciK* , (combined) 



18 6 ± 2.4 0.20 0.097 ± 0.040 ± 0.009 

16 1±1 0.14 < 0.27 
15 1.2 ±1.5 0.13 < 0.21 
15 0.033 < 0.67 

17 0.11 < 0.30 
< 0.21 



a This efficiency does not include the ip,ip ,x c i,K ,K* or Kg branching ratios. The two sets of values 
given for the ip channels correspond to the two ip decay modes ip and ip ipir + ir~ . 
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TABLE X. Upper limits (90% C.L) on branching fractions for color suppressed B decays. 



T) prav 1\/T nnp 
i / \ \ <i y _lv_i_u u_c 


KvPTl t C! 

J_J V Cll Li O 




TI T, (%) at 90% C T, 


RO no o 


^ 90 7 
\ zu. 1 


n 39 n 1 r n 1 s 




RO nO 


/ion 


PI 91 (1 (18 (119 


\ vJ.vJOO 


ro > no 


/ r, 

\ .0 


n 31 n 1 1 n 1 


\ u.uuo 


2?° -+ I} ?/ 


< 3.5 


0.18,0.08,0.11 


< 0.086 


2?° -+ 


< 12.7 


0.16,0.07,0.09 


< 0.063 


5° _D* 7T° 


< 11.0 


0.13,0.07,0.07 


< 0.097 




< 8.1 


0.09,0.04,0.04 


< 0.117 




< 2.3 


0.11,0.05,0.06 


< 0.069 




< 2.3 


0.07,0.03,0.03 


< 0.27 


50 _, £)*0^ 


< 9.0 


0.06,0.03,0.03 


< 0.21 


e The efficiencies 


for the D° -> K~tt + , D° - 


-> K-tt+tt , and D° - 


-> K~ir + ir~'ir + modes are given. 



These efficiencies do not include D, D* r], r\ and u branching ratios. 



TABLE XL Upper limits on ratios of branching ratios for color suppressed to normalization 
modes. 



Ratio of Branching 


Ratios 


U.L. (90% C.L.) 


B(B° - 


■* D Tr°)/B(B- 


D°ir-) 


< 0.09 


B(B° - 


-> D°p°)/B(B- 


- D°p-) 


< 0.05 


B(B° - 


-* D°r])/B(B~ 


-+ D°TT-) 


< 0.12 


B(B° - 


-> D T]')/B(B- 


-+ D°TT-) 


< 0.16 


B(B° - 


-> D°u)/B(B- 


- D°p-) 


< 0.05 


B(B° ^ 


_D*°7T )/ B(B~ 


D*°ir~) 


< 0.20 


B(B° - 


D*°p )/B(B- 


-+ D* p~) 


< 0.07 


B(B° - 


» D*°r])/B(B~ 


-+ D*°TT~) 


< 0.14 


B(B° - 


► D*°r]')/B(B~ 


-+ D*°tt-) 


< 0.54 


B(B° - 


> D*°u)/B(B- 


-+ D*°p~) 


< 0.09 



TABLE XII. B Masses from individual modes (not corrected for initial state radiation). 



Mode 


B~ Modes 
Mass (MeV) 


Events 


Mode 


i? Modes 
Mass (MeV) 


Events 


D*°ir~ 


5279.7±0.4 


73 


D*+ir- 


5280.1±0.4 


73 


D*°p~ 


5280.2±0.4 


89 


D*+p~ 


5280.5±0.4 


79 




5279.8±0.4 


44 


ij)K*° 


5280.4±0.5 


29 


D°TT- 


5279.9±0.3 


76 


D+ir- 


5280.4±0.3 


80 


D°p~ 


5279.7±0.4 


80 


D+p- 


5280.3±0.4 


79 


All 


5279.9±0.2 


362 


All 


5280.3±0.2 


340 
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TABLE XIII. Contributions to the systematic error in the B° — B mass difference. 

Event sampie 0.15 MeV 

Background shape 0.12 MeV 

7 energy caiibration < 0.03 MeV 

Width of B mass peak < 0.02 MeV 

Track momentum scaie < 0.01 MeV 

total 0.19 MeV 



TABLE XIV. Measurements of the B°-B~ Mass difference. 



Experiment M(B°)-M(B-) (MeV) 

CLEO 87 [2] 2.0 ± 1.1 ±0.3 

ARGUS [6] -0.9 ±1.2 ±0.5 

CLEO 92 [3] -0.4 ±0.6 ±0.5 

CLEO 93 (this result) 0.41 ± 0.25 ± 0.19 

Average 0.4 ± 0.3 



TABLE XV. Ingredients for Factorization Tests. 



|ci| 1.12 ±0.10 

U 131.74 ± 0.15 MeV 

f p 215 ± 4 MeV 

f ai 205 ± 16 MeV 

V ud 0.975 ±0.001 

j^(B -> D*l v)\ q i =ml {WSB) 0.0023 GeV" 2 

j^{B -+ D*l v)\ q2=m 2 (I SGW) 0.0020 GeV" 2 

9 j^(B -+ D*l v)\ q J ml (KS) 0.0024 GeV" 2 

j^(B D*l v)\ q 2 =m 2{WSB) 0.0025 GeV" 2 

j^(B D*l v)\ q 2 = J p {ISGW) 0.0024 GeV" 2 

^{B -+ D*l v)\ q 2= m 2 p {K S) 0.0027 GeV" 2 



dq 

j^(B -> D*l v)\ q 2 =m 2^ (WSB) 0.0032 GeV 

dt 

dq 



OL( n v n*7 ,,M „ „ ("wen"! n rwwo n n u-2 

j^(B D*l v)\ q 2 =m 2 a \lSGW) 0.0030 GeV" 2 

j^(B D*l v)\ q 2 = ^2 ai (KS) 0.0033 GeV" 2 



TABLE XVI. Comparison of R exp an( i Rth< 







R exp (GeV 2 ) 


Rtheor (GeV 2 ) 


B° - 


-> D*+ir- 


1.1 ± 0.1 ± 0.2 


1.2 ± 0.2 


B° - 




3.0 ± 0.4 ± 0.6 


3.3 ± 0.5 


B° - 


-> D*+a^ 


4.0 ±0.6 ±0.5 


3.0 ± 0.5 
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TABLE XVII. Ratios of B decay widths. 



Exp. Factorization RI Model BSW Model 

B(B° D*+p-)/B(B° D*+tt-) 2.9 ±0.5 ±0.5 2.9 ± 0.05 2.2 - 2.3 2^8 
B(B° -r D* + a-)IB(B -r D* + tt-) 5.0 ±1.0 ±0.6 3.4 ± 0.3 2.0-2.1 3.4 



TABLE XVIII. Ratios of color suppressed to external spectator branching ratios. 

B(D° -+ K°ir )/B(D° -+ I(-ir+) 0.57 ± 0.13 

B(D° K* Tr )/B(D° -> K*~tt + ) 0.47 ± 0.23 

B(D° -+ Tr Tr )/B(D° -+ 0.77 ±0.25 

B(D+ -> K*°K+)/B(D s -> <M + ) 0.95 ± 0.10 

B(D+ -> K°K+)/B(D s -> <f>ir+) 1.01 ± 0.16 



TABLE XIX. Branching fractions of color suppressed B decays and comparisons with models. 



Decay Mode 


U. L. (%) 


BSW (%) 


B (BSW) 


B° - 




< 0.048 


0.012 


0.20a^(/ D /220MeV) 2 


B° - 


-> D°p° 


< 0.055 


0.008 


0.14a 2 (/ D /220MeV) 2 


B° - 


-> D°T] 


< 0.068 


0.006 


0.11a2(/ D /220MeV) 2 


B° - 


-> D° V ' 


< 0.086 


0.002 


0.03a 2 (/ D /220MeV) 2 


B° - 


-> D°uj 


< 0.063 


0.008 


0.14a 2 ,(/ D /220MeV) 2 


B° - 




< 0.097 


0.012 


0.21a^(/ D */220MeV) 2 


B° - 


-> D*°p° 


< 0.117 


0.013 


0.22a 2 (/ D */220MeV) 2 


B° - 


-> D*°T] 


< 0.069 


0.007 


0.12a 2 (/ D */220MeV) 2 


B° - 


-> D*°T]' 


< 0.27 


0.002 


0.03a 2 (/ D */220MeV) 2 


B° - 


■+ D*°cj 


< 0.21 


0.013 


0.22a^(/ D */220MeV) 2 



RI model(%) 



0.0013 - 0.0018 
0.00044 



0.0013 
0.0013 



0.0018 
0.0018 
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TABLE XX. Branching ratios in terms of BSW parameters a\, a 2 



Mode 






B % 


B° -+ 


D+ir- 




0.264a 2 


B° -+ 


D+p- 




0.621a 2 


B° 


D* + ir~ 




0.254a 2 


B° -+ 


D*+p~ 




0.702a 2 


B~ -> 


■ D°TT- 


0.265K + 


1.230a 2 (jb/220)] 2 


B~ -> 


■ D°p- 


0.622[ai + 


0.662a 2 (/ D /220)] 2 


B~ -> 


- D*°ir~ 


0.255K + 


1.292a 2 (/ D ./220)] 2 


B~ -> 


- D*°p~ 


0.703[a 2 + 0.635a 2 (/ D 


./220) 2 + 1.487aia 2 (/d*/220)] 


B~ -> 


■ ipK~ 




1.819a| 


B~ -> 


- ipK*~ 




2.932a^ 


B° -+ 


ij)K° 




1.817a 2 


B° -+ 






2.927a 2 



TABLE XXI. Ratios of normalization modes to 


determine the sign of a 2 /ai. 


The magnitude of 


a 2 /ai is 


the value in the BSW 


model which agrees 


with our result from B if) 


modes. 


Ratio 


a 2 /a 1 = -0.24 


a 2 /a 1 = 0.24 


CLEO II 


RI model 


Ri 


0.50 


1.68 


1.89 ±0.26 ±0.32 


1.20 - 1.28 


R2 


0.71 


1.34 


1.67 ± 0.27 ± 0.30 


1.09 - 1.12 


Rs 


0.48 


1.72 


2.00 ±0.37 ±0.28 


1.19 - 1.27 


i?4 


0.41 


1.85 


2.27 ± 0.41 ± 0.41 


1.10 - 1.36 
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XIV. APPENDIX 



In this appendix, we provide the product of the B and charm branching fractions for the 
decay modes measured in this paper so that the results can be easily renormalized when the 
intermediate branching fractions for D°, D + , D* + , D*° and ip, if) , Xci mesons are known more 
precisely. The results are given in Tables XXII-XXVIII. 



TABLE XXII. Product Branching Fractions (%) for B~ D°(mr)- Modes 



B~ Mode 


D Mode 


#of 


e 


B(B~ -+ D°(nir)-)x 






events 




B(D° -+ K-[mr]) 




K~ir + 


76.3±9.1 


0.433 


0.0189 ± 0.0022 ± 0.0013 


D°TT- 


K-ir+ir 


134±15 


0.193 


0.0746 ± 0.0082 ± 0.0065 




K~ir + ir + 'ir~ 


94± 11 


0.222 


0.0455 ± 0.0055 ± 0.0049 




K~ir + 


80 ± 9 


0.155 


0.0524 ± 0.0067 ± 0.0044 


D°p~ 


K-ir+ir° 


42 ± 9 


0.036 


0.1254 ± 0.0282 ± 0.0150 




K~ir + ir + 'ir~ 


90.4±12.1 


0.079 


0.1223 ± 0.0164 ± 0.0142 



TABLE XXIII. Product Branching Fractions (%) for B° ^"^(raTr)- Modes 



B° Mode 


D Mode 


#of 

events 


e 


B(B° -+ D+(nir)-)x 
B(D+ -> K-ir + ir + ) 


D+ir- 


K~ir + ir + 


80.6±9.8 


0.32 


0.0265 ± 0.0032 ± 0.0023 


D+p- 




78.9±10.7 


0.12 


0.0704 ± 0.0096 ± 0.0070 



TABLE XXIV. Product Branching Fractions (%) for B~ D*°(mr) Modes 



B~ Mode 


D° 


Mode 


#of 

events 


e 


B(B~ -+ D*°(nir)-)x 
B(D*° -+ D°ir°) x B(D° -+ K-[mr]) 




K~ 


■7T + 


13.3±3.8 


0.16 


0.0090 ± 0.0026 ± 0.0009 


D*°ir~ 


K~ 


■7T + 7T 


37.7±6.9 


0.08 


0.0488 ± 0.0089 ± 0.0063 




K~ 


'7r + 7r + 7r~ 


20.0±4.9 


0.08 


0.0267 ± 0.0065 ± 0.0033 




K~ 


■7T + 


25.7±5.4 


0.064 


0.0432 ± 0.0090 ± 0.0058 


D*°p~ 


K~ 


■7T + 7T 


43.8±7.8 


0.027 


0.1722 ± 0.0305 ± 0.0300 




K~ 


'7r + 7r + 7r~ 


16.9±4.6 


0.030 


0.0608 ± 0.0176 ± 0.0095 




K~ 


■7T + 


5.5±2.9 


0.048 


0.0124 ± 0.0065 ± 0.0020 


D*°ir-ir-ir +a 


K~ 


■7T + 7T 


27.7±7.2 


0.022 


0.1316 ± 0.0343 ± 0.0237 




K~ 


'7r + 7r + 7r~ 


15.0±4.5 


0.025 


0.0632 ± 0.0187 ± 0.0118 



a The three pion mass is required to be between 1.0 GeV and 1.6 GeV consistent with an (i\ meson. 
(If this channel is dominated by aj", the branching ratio for _D*°a 1 " is twice that for _D*°7r _ 7r~7r + .) 
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TABLE XXV. Product Branching Fractions (%) for B° D*+(mr)- Modes 



B° Mode 


D° Mode 


#of 


e a 


B(B° -+ D*+(nir)-)x 






events 




B(D*+ -> D°ir+) x B(D° -> K-[mr]) 




K~ir + 


19.4±4.5 


0.35 


0.0058 ± 0.0013 ± 0.0008 




A 7T 7T 


91 0_L£ /I 

oi.y±o.4 


n 1 a 
U.14 


U.UZ4o ± U.UU4y ± U.UUoO 




A"~7r + 7r + 7r~ 


20.5±5.2. 


0.15 


0.0146 ± 0.0033 ± 0.0025 




A"~7T + 


21.9±5.2 


0.12 


0.0188 ± 0.0044 ± 0.0034 






39.8±7.2 


0.048 


0.0892 ± 0.0162 ± 0.0177 




li' _ 7r + 7r + 7r~ 


14.6±4.6 


0.054 


0.0286 ± 0.0091 ± 0.0059 




A"~7T + 


13.5±3.9 


0.096 


0.0151 ± 0.0044 ± 0.0024 


J D*+7r-7r-7r +a 


K-ir+ir 


21.7±5.9 


0.043 


0.0545 ± 0.0147 ± 0.0091 




li' _ 7r + 7r + 7r~ 


13.9±4.4 


0.042 


0.0348 ± 0.0101 ± 0.0069 


a The three pion mass is required to 


be between 1.0 GeV 


and 1.6 GeV consistent with an a\ meson. 


(If this channel dominated by a 1 , the branching 


ratio for D* + a 1 is twice that for D* + ir it 7r + .) 


TABLE XXVI. 


Product Branching Fractions for B 


— ► ip Modes and 90% Confidence Level 


Upper Limits (%). 










B Mode 


# of events 


e 




B(B -+ ipK^) x B(ip -+ /+/-) a 




58.7 ± 7.9 


0.47 




0.0131 ±0.0017 ±0.0011 


B° 


10.0 ±3.2 


0.34 




0.0088 ± 0.0028 ± 0.0009 




29.0 ±5.4 


0.23 




0.0200 ± 0.0037 ± 0.0021 


5" VA'*" 


12.6 ±3.6 






0.0210 ± 0.0061 ± 0.0026 



a The product branching fraction has been corrected for the K°, K* or Kg branching ratios but not 
for the ip branching fractions. 
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TABLE XXVII. Product Branching Fractions for B ip Modes and 90% Confidence Level 
Upper Limits (%). 



B Mode 

B~ ^'K-if}' -> /+/- 
B~ — ► ip K~,ip — ► '07r + 7r _ 

B° ij;'K ,ij}' -> ^tt+ti-- 
5° tl}'K*°,4>' -> /+/- 
5° tl}'K*°,4>' -> ^tt+ti- 
5" i,'K*-,i) -> /+/- 
_B~ — ► -0 K*~,ip — ► -(/>7r + 7r _ 



0.0011 ± 0.0006 ± 0.0001 
0.0173 ± 0.0100 ± 0.0023 

< 0.0025 

< 0.0200 

< 0.0051 

< 0.0210 

< 0.0065 

< 0.0600 



a The product branching fraction has been corrected for the K°, K* or branching ratios but 
not for the ip and ip branching fractions. We give B(B ip K^) X B(ip l + l~) or B(B 

iP'l&y) X B(lp' -+ IpTT + TT-) X B(lp -+ / + /-). 



TABLE XXVIII. 


Product Branching 


; Fractions for B - 


-> Xd Modes and 90% Confidence Level 


Upper Limits (%). 








5 Mode 


# of events 


e 


B(B -+ XciKM) x B(xci) a 




6 ± 2.4 


0.20 


0.0031 ±0.0013 ± 0.0003 




1 ± 1 


0.14 


< 0.0087 


B° -+ XciK* 


1.2 ± 1.5 


0.13 


< 0.0066 


B~ -+ XciK*- 







< 0.0066 



a For the modes with Xd mesons, we report the product B(B [x c i]K^)xB(ip l + l~)><Xci - ► ~fp 
branching fraction in the product. The product branching fraction has been corrected for the K°, 
K* or Kg branching ratios but not for the ip branching fractions. 
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